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We present the so-called Liouville-Neumann (LN) approach 
to nonequilibrium quantum fields. The LN approach unifies 
the functional Schrodinger equation and the LN equation for 
time-independent or time-dependent quantum systems and 
for equilibrium or nonequilibrium quantum systems. The LN 
approach is nonperturbative in that at the lowest order of 
coupling constant it gives the same results as those of the 
Gaussian effective potential at the zero and finite tempera- 
ture in a Minkowski spacetime. We study a self-interacting 
quantum field in an expanding Friedmann-Robertson- Walker 
Universe. By studying a toy model of anharmonic oscillator 
and finding the underlying algebraic structure we propose a 
scheme to go beyond the Gaussian approximation. 



I. INTRODUCTION 

Recently there have been much interests in the cos- 
mological scenarios. These models have been developed 
to describe the evolution of the Universe from the early 
stage to the present epoch ||]. Quantum effects or high 
temperature thermal effects of matter fields have dom- 
inated in the early Universe. But in the cosmological 
scenarios the Universe has been evolving and expanding. 
Quantum fields in such an expanding background space- 
time could not be described properly as thermal equilib- 
rium. They are characterized by nonequilibrium quantum 
processes. Quantum nature of the matter fields predom- 
inated as it dated back into the early stage of the Uni- 
verse. Several field theoretical methods have been intro- 
duced to describe properly the nonequilibrium quantum 
process. Typical and frequently used methods are closed 
time path and thermal field dynamics [0 . Closed time 
paths integrate over contours of imaginary (or complex) 
time and real time. Thermal properties of quantum field 
are exhibited through the contours of imaginary time. 
Thermal field dynamics doubles the degrees of freedom of 
the system and thermal states are condensation of parti- 
cles of the system and its fictitious counterpart. Despite 
wide applications closed time paths are difficult to do 
integrals for self-interacting quantum fields and to incor- 
porate various initial quantum conditions, ft would be 
more than just an interest to have another method, par- 
ticularly a canonical method, for nonequilibrium quan- 
tum fields. 



In this paper we present another field theoretical 
method, the so-called Liouville-Neumann approach, to 
describe the nonequilibrium quantum processes in cos- 
mology [|]]^]. It is a nonperturbative canonical method 
that unifies the Liouville-Neumann (LN) equation and 
the functional Schrodinger equation H. It is nonpertur- 
bative in that at the lowest order of coupling constant 
the wave functionals take into account nonperturbative 
quantum effects and are equivalent to an infinite summa- 
tion of daisy and superdaisy diagrams or the Gaussian 
effective potential. In particular, it is the LN equation 
that describes the nonequilibrium quantum process. The 
underlying assumption of the LN approach is that all 
nonequilibrium processes are consequence of microscopic 
processes which are well described by quantum theory. 
An introduction of canonical method will enlarge our 
understanding of the nonequilibrium quantum processes 
that have mostly been dealt with the closed time path in- 
tegral methods but are still not completely understood. 
Furthermore, the LN approach has several advantageous 
points compared with the path integral methods. First 
of all, it is truly canonical. Many useful concepts and 
ideas from quantum mechanics can straightforwardly be 
applied to quantum field theory by treating a quantum 
field as a system with infinite degrees of freedom. Sec- 
ondly, the Hilbert space of the exact or approximate wave 
functionals is constructed explicitly. Therefore, it is rel- 
atively easy to incorporate various initial quantum con- 
ditions and to evolve them. Moreover, the LN approach 
provides a scheme to go beyond the Gaussian approxi- 
mation. Thirdly, the LN approach can readily be unified 
with thermal field dynamics due to the canonical nature. 
Besides these there are also some other useful points. 

The organization of this paper is as follows. In 
Sec. II we introduce the LN equation to solve time- 
dependent quantum systems and unify it with the func- 
tional Schrodinger equation to describe quantum fields. 
In Sec. Ill the LN approach is applied to a mas- 
sive scalar field in an expanding Friedmann-Robertson- 
Walker (FRW) Universe and the quantum states or 
Wightmann functions are found for various initial con- 
ditions. In Sec. IV we apply the LN approach to a 
self-interacting quantum field in the FRW Universe and 
find various quantum states or wave functionals. A toy 
model of anharmonic oscillator is studied to go beyond 
the Gaussian approximation. Comparison of the LN ap- 
proach is made with the Gaussian effective potential. In 
Sec. V the LN approach is further exploited beyond the 
Gaussian approximation and reveals the algebraic struc- 



f 



ture of the anharmonic oscillator. 



B. Functional Schrodinger-Picture 



II. LIOUVILLE-NEUMANN APPROACH 

As mentioned in the Introduction there are several mo- 
tivations to study the LN approach to quantum field 
theory. The LN approach can be used in a unified way 
to both time-independent and time-dependent quantum 
systems, and to both equilibrium and nonequilibrium sys- 
tems. It is a canonical method that deals with the non- 
perturbative quantum effects. It is canonical because it 
provides a method to solve the functional Schrodinger 
equation. As a field theoretical method we wish to unify 
the LN equation and the functional Schrodinger-picture 
to make the so-called LN approach. 

A. Liouville-Neumann equation 

The LN approach to quantum mechanics has been 
introduced by Lewis and Riesenfeld |?J to solve time- 
dependent quantum systems, especially time-dependent 
harmonic oscillators. It is based on the observation that 
the Schrodinger equation (h = 1) 



(i) 



for time-dependent as well as time-independent systems 
can be solved in terms of the following operator 



d 



+ [6,H}-i^-6 = 0, 
V 



0\Vn(t))=Vn(t)\Vn(t)). 

The exact quantum states are indeed given by 



(2) 



/ / dt(r] n (t)\i^- t - JIM*)) //„!/);• 



(3) 



The functional Schrodinger-picture approach to quan- 
tum field theory is a canonical method in which the wave 
functionals of a field as quantum states are determined 
by the functional Schrodinger equation B 



i-*(ct>,t) = H( ( t>,-i— )*(<M), 



(5) 



where <p represents the quantum field. The Hamiltonian 
describes an infinite dimensional system since the field 
has infinite degrees of freedom. The wave functionals 
^>((p,t) are specified by assigning c-numbers to the field 
</>(x) at a fixed time. The wave functionals constitute the 
Hilbert space equipped with the inner product 



<*i|tf 2 > = / X>[fl*;(0)* a (0). 



(6) 



Operators act on the wave functional of the Hilbert space 



0(0, tt)|*(0, tt)) 0(0, -i—)^{cj>). 



(7) 



The Gaussian wave functionals are a generalization of the 
Gaussian states of a harmonic oscillator: 



* o (0) = iVexp 



(<Mx)-</>(x,i)) 

x.y 

-iE(x,y,i))(0(y)-0(y,i)) 

(8) 



.4G(x,y,t) 

+1 { 7f(x,t)(0(x)-0(x,t)) 



where 



<*o|0(x)|*o>=0(x,t), 

(*0|7T(X)|*0> =7f(x,t), 

(*o|0(x)0(y)|*o) = 0(x, W(y, t) + G(x, y, t), 
(*o|7r(x)7r(y)|*o) = 7f(x, t)7f (y, t) + E(x, y, t). (9) 



In particular, when O satisfies the LN equation 

ild + [d,A\ = o, 



(4) 



its eigenvalues are time-independent, rj n — 0. Such an 
operator is called the Lewis- Riesenfeld invariant operator 
or generalized invariant operator or LN operator. 

The aim of this paper is to develop this LN approach 
to quantum mechanical systems as a field theoretical 
method in the functional Schrodinger-picture. 



III. FREE SCALAR FIELD IN THE FRW 
UNIVERSE 

As an application of the LN approach, we consider a 
free massive scalar field in a spatially flat FRW Universe 
with the metric 



ds 2 = -dt 2 + R 2 (t)dx 2 . 
The scalar field has the Lagrangian 

i? 3 n 2 (V0) 2 



(10) 



L 



d J x- 



R 2 



(m 2 + £ (4) R)^ , (11) 
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where £ = 0, g for the minimal and conformal couplings, 
respectively, and ( 4 -*R is the scalar curvature. The Hamil- 
tonian has the form 



H = I d 3 x 



^ + |W) 2 + f (™ 2 +^>R)^ 

(12) 

We quantize the scalar field according to the functional 
Schodinger-picture 



(13) 



Knowing the Fock space of exact quantum states for 
a time-dependent harmonic oscillator, we decompose the 
Hamiltonian into a sum of harmonic oscillators. For this 
purpose, we decompose the field into Fourier-modes 



(x,t) = 



1 



(2tt) 3 / 2 

and refine the field modes as 



d 3 kcf) k (t)e 



(14) 



4 +) = ^(^k + 0-k), 



(15) 



Space integrations can easily be done for the field and 
momentum: 

d 3 X0 2 (x,i) = f d 3 fc0 k 0_ k = ^0 2 , 

f d 3 X7T 2 ( X ,t) = f d 3 fc7r k 7T_ k EE ^TT 2 ,, (16) 

ot 

where a denotes (±, k). Then the Hamiltonian takes the 
form 



R 3 u 2 



L2i? 3 



where 



^W = ™ 2 + ^+e (4) R- 



(17) 



(18) 



That is, the Hamiltonian ( |l7| ) is a countably infinite 
sum of time-dependent harmonic oscillators. The func- 
tional Schrodinger equation ( |l3| ) becomes the ordinary 
Schrodinger equation for the infinite system 

a 

The wave functional is a product of each wave function 
*(0)=JJ*(&,). (20) 



In the end the right hand side will be expressed in terms 
of the field </> by taking an appropriate inverse Fourier- 
transform. 

In the LN approach, one looks for the following first 
order operators 



&1 = -i(ip a (t)lT a - <p a (t)4> a ). 

And then one imposes the LN equations 
d 

i-Q^a a + [a a ,H a ] = 0, 
d 

i-Q t ai+[al,H a ] =0. 



(21) 



(22) 



Equations (|22j) are satisfied when each mode satisfies the 
classical equation of motion 



R 

'fait) + 3 — tp a (t) + LU 2 a Lp a (t) = 0. 

H 



(23) 



One further requires a a and a) a to form the basis for the 
Fock space: 



(24) 



These commutation relations lead to the boundary con- 
ditions on ip a 



(25) 



Then the Fock space of the number states of each mode 
is constructed by 



Oa|0 Q ,t> = 0, 

1 f \ n <* 

\n a ,t) = . — - (a a ) \0 a ,t}. 



(26) 



The vacuum state of the field is the product of vacuum 
states for each mode 



|0,t)=ni°a.*>- 



(27) 



A few comments are in order. The vacuum state ( |27| ) has 
all the symmetries of the Lagrangian. Knowing the Fock 
space of exact quantum states, it is easy to incorporate 
various initial quantum conditions. 

Usually initial conditions are prepared at a particular 
time to- The vacuum and the number states are denoted 
by |0 Q ,to) and \n a ,ta). First, with respect to the initial 
vacuum state, one finds the Wightmann function 



<0(x,t) ^(X',t'))vac. 



(2tt) 



— y 



a i(ki-x— k 2 -x') 



vac. 



ai .CK2 

d 3 k 



(2tt) e 



«k- (x— x') * 



vUt)<Pk>(t')- 



(28) 
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Second, for an initial thermal state 



A. (0 + l)-Dimensional Toy Model 



Pa (*o) 



1 



-PH a (t ) 



the Wightmann function is found 



(29) 



d 3 k 
(2^)3 



Jk-(x-x') 



1 



x(^W^(0 + ^(iM-(t')). (30) 

Finally, let us consider a coherent state described by the 
density operator 



p a ji(to) — 



1 



Trp aJI (t 



_ 7. 



where e a 
thermal state displaced by 



< \p - fi(u a (t Q )a) a a a 

+ 7q«1 +7 q «q +e«) i (31) 
Equation ( |3l| ) indeed describes a 



r< + 7 



D a = e -»('o)^ 1 "o(* r 
The Wightmann function is found 

(0(x,^(x',t'))coh. = (0(x,i)0(x',t')}thcrm. 



(32) 



(2tt) 3 



¥>k,c(*)<Pk' c (*') 5 



where 



<Pk, c = — (7k<P k + 7k<Pk), 



(33) 



(34) 



is a classical field. 



IV. SELF-INTERACTING SCALAR FIELDS 

In this section we apply the LN approach to a self- 
interacting quantum field. There have been intro- 
duced different methods to investigate the nonperturba- 
tive quantum effects. The most typical nonperturbative 
method is the Gaussian effective potential (variational 
Gaussian approximation) approach But this method 
has been applied to time-independent quantum fields. 
The Gaussian wave functionals have also been considered 
in Ref . |j] . The problem of solving the covariant kernel is 
not completed. Compared with these methods, the LN 
approach provides us with the explicit wave functionals 
in terms of the field equations satisfied by each mode. 
The LN approach has recently been applied to quantum 



fields in curved spacetime 



ntiy b 



To illustrate how the LN approach works for a self- 
interacting quantum field, first we consider a time- 
dependent anharmonic oscillator. More than just a quan- 
tum mechanical system, the anharmonic oscillator has 
been used as an important toy model to test various field 
theoretical methods. 

Let us consider the anharmonic oscillator 



H = 



3 2 , ^{t) 
2 + — q 



— q - 



(35) 



The mass has been rescaled for convenience. As in the 
case of harmonic oscillators in Sec. Ill, at the lowest 
order of the coupling constant one optimizes the Hilbcrt 
space by a Fock space of the annihilation and creation 
operators 



a(t)=i{y*(t)p-y*(t)q), 
tf(t) = -i(y(t)p-y{t)q). 



(36) 



We require these operators to minimize the energy or 
the Dirac action in the time-independent case or time- 
dependent case, respectively. In the time-dependent case 
we can not use the minimization of the energy, because 
the energy is not conserved. In the time-independent 
case, these two principles of minimization give the iden- 
tical result. The resulting quantum state is the Gaussian 
wave function. In the time-dependent case, it has been 
found that the LN approach gives the same result as that 
of the Dirac action principle at the lowest order of the 
coupling constant. However, the LN approach has sev- 
eral advantageous points. First, it is an algebraic method 
for the anharmonic oscillator problem. The underlying 
group structure proves very convenient and useful. The 
underlying structure of the anharmonic oscillator will be 
studied in detail in the next section. Second, the quan- 
tum states can be found nonpcrturbatively to any desired 
order of the coupling constant. 

Then the anharmonic oscillator has the oscillator rep- 
resentation 



H = H G + XH' 



(37) 



where 
H G = 



y*y + uj 2 y*y + 3\(y*y) 

1 

~2 

2 , , .2. .2 , o \ f„.*„.\„2 



at " + 2 1 ) ~T (y * y)2 



y* +uj y* +3\{y*y)y* 



,t2 



4 ^ \k 

k=0 



y *(4-fc) y fc a t(4-fc) a fc. 



(38) 
(39) 



At the lowest order of the coupling constant, we use the 
approximate LN equations 



4 



= 0, 
= 0. 



Equations ( p0| ) lead to the equation for y 

y(t)+u J 2 (t)y + 3\(y*(t)y(t))y(t) = 0. 



(40) 



(41) 



The state annihilated by a becomes an approximate 
ground state 

a|0) [0] - 0, (42) 
and in the coordinate representation is given by 

' ' " ' (43) 



v 1/4 


r- v* 2] 


j exp 


i q 

- 2y* i 



Another advantageous point of the LN approach is that 
the excited states can be found easily 



1 



*[o],n(g.*) = -H(a t ) *[o],o(?,*) 



(44) 



A few comments are in order. Though approximate, 
these quantum states take into account nonperturbative 
quantum effects already. In field theory, the Gaussian 
state (^) is equivalent to a summation of daisy and su- 
perdaisy diagrams. For a time-independent harmonic os- 
cillator the Gaussian state (43) is also known to account 
for the energy within a few percent even for the strong 
coupling limit. 

We now wish to find the operators that satisfy the LN 
equations to any desired order. By treating H' pertur- 
batively in the full LN equations 

8 



ig- t A+[A,H]=0, 



(45) 



one expands these operators in a series of A 

oo 
oo 

^ t =« t + E A "A«)- 



(46) 



In the Fock space representation, B/. 
panded in a and a^: 



B, 



J2 b 



(n) a a ■ 



By comparing the same powers of A in Eq. 
the recursive relations 

d - * 

— B {n) = i [B [n) ,H G ]+i [B(„-i) , H'} . 



are again ex- 

(47) 
3), we find 



(48) 



The recursive relations can be solved iteratively 
B(i)=< f [a,H'], 



B, 



>(n)=i n J ■■ J[[---[[a,H'],H'],...,],H']. (49) 
For instance the first order corrections are found 
B (1) =4(i y y* 4 )at3 + 3 (z J y* 3 y)a^a 

+2( i | ? /V)« t « 2 + (*/yV)« 3 . (50) 

An improved ground state is obtained from 

Ai]l°>[i] = °> ( 51 ) 

where 

A m =a + XB {1) . (52) 

This procedure can be repeated to solve Eq. ((45|). The 
quantum states at any order of the coupling constant are 
obtained therefrom. 

A comment is in order. In the language of field theory, 
the wave functional determined by a and a) corresponds 
to a free propagator with a renormalized mass and that 
determined by A and ^4^ m (^6|) corresponds to a 
summation of all loop diagrams. In this sense Eq. ( |49| ) 
can be interpreted as the Feynman rules on the closed 
time path for the time-dependent system. 

B. Self-Interacting Scalar Field 

We now turn to a self-interacting quantum field. A 
quantum field is equivalent to a system of infinitely mode- 
decomposed harmonic or anharmonic oscillators. Thus 
the results from the anharmonic oscillator system are ex- 
pected to be useful in studying the self-interaction quan- 
tum field. Though the formalism (LN approach) to be 
developed in this paper can be applied to any quantum 
field theory, either renormalizable or nonrenormalizable, 
we shall confine our attention to the (/) 4 -theory. The rea- 
son is that the results from our formalism can readily be 
compared with those from other methods. 

The 4 -theory in the FRW Universe has the Hamilto- 
nian density 



H 



R 



/m 2 



A 



(53) 



Before the Fourier-decomposition, the field is divided into 
a classical background field and quantum fluctuations 



4> — 4> c + 0/- 



(54) 
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The classical background field depends only on the co- 
moving time and equals to the space average of the field 
(zero-mode). Quantum mechanically it is a coherent 
state of the field 



{d>) = d>M 



(55) 



We assume the fluctuations to have symmetric quantum 
states 



(h) = o. 



(56) 



Then the Hamiltonian density is the sum of classical 
background field part and quantum fluctuations part: 



H = 



+R 3 ( m2 3A 



a 



"m 

Y 



+ 



XR 3 



R? 



R 



3 W 



(57) 



The expectation values of the terms in the last square 
bracket vanish with respect to the symmetric quantum 
states of fluctuations. So these terms will not considered 
any more. 

First, we study the classical background field. Only 
those terms whose expectation values do not vanish with 
respect to the symmetric quantum states of fluctuations 
contribute to the Hamiltonian density for the classical 
background field 



H B = 



2R 3 



„o/m 3A , - 2 



(58) 



The classical background field is frequently treated as a 
classical system. Then the equation of motion is 

4>c + 3^ c + (m 2 + 3\{4>}))(t> c + Xf c = 0. (59) 

The infinite contribution from the quantum fluctuations 
should be regulated in a suitable manner. But more 
proper treatment is to quantize even the classical back- 
ground field. Therefore, the background field is nothing 
but an anharmonic oscillator with the mass shifted to 
to 2 + 3A(0 2 ). The classical background field plays the 
role of an inflaton in the inflationary scenarios. The im- 
portance of such a quantum background field has been 
shown in Rcf . [ . The annihilation and creation opera- 
tors of the Fock space are 

d c = i(ip* c TT c - R 3 ip c *(j) c ) - i(ip* c n c - R 3 ip c *4> c ), 

a\ = -i(ip c (t)TT c - R 3 (p c 4> c ) + i(ip c n c - R 3 ip c <j) c ). (60) 

The d c and dj are chosen so that the vacuum state gives 
the expectation values 



R 3 



(61) 



Requiring the LN equations for d c and dj, it is shown 
that ip c satisfies the equation 

R 

</3 c + 3-y> c + (m 2 + 3\{4> 2 f ))ip c + \(ip* c (p c )ip c = 0. (62) 

In the abo ye p rocedure, <p c is required to satisfy the same 
equation ( p2|) as ip c - So the terms proportional to the 
identity operator in Eq. ( ]60| ) become the Wronskian of 
the equation of motion, and constant. This makes our 
procedure self-consistent. In the minimal uncertainty 
proposal for the inflation jl(J, a particular choice of co- 
herent state was assumed 



0. 



(63) 



and the quantum states of the inflaton were again the 
symmetric Gaussian states. 

Next, we consider the quantum fluctuations. The 
terms in the second square bracket of Eq. ( |57j) give rise 
to the Hamiltonian density for quantum fluctuations 



Hi 



2R 3 



+ | ( V0 /) 2 + ^(^ + ^ 2 )0 2 + 



f 



(64) 



Upon decomposing into Fourier-modes, we keep only 
those terms with nonvanishing contribution 



(65) 



We express the fluctuations and their momenta in terms 
of the annihilation and creation operators 



4> a = ¥> a a a + y*d^, 

Ti a = R 3 (tp a a a + <Paa£) ■ 



(66) 



As for the anharmonic oscillator, we find the oscillator 
representation of the Hamiltonian for fluctuations 



H f = H,, G + XH' 



(67) 



where -ff/,G and H'j are the quadratic and non-quadratic 
parts. At the lowest order of the coupling constant we 
are interested in Hf^c- 



Hf.c = R 3 ^2 

a 

+ 3A 



((2tt) 3 



(27T) 

3Ai? 3 1 
4 (2tt) 



<PpVp)<P a <Pa 



3A0 2 )(^ Q 

(ala a + - 
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R 3 ^ 



( k 2 



I R 2 



*t2 



R 3 x 



+ 3A ((^E^K' 



/3 



(68) 



The equations of motion for the variables <p a 
R /ok 2 



+3A ((^3X]^)^ = ' 



(27T) 

are determined from the LN equations 
d 



(69) 



(70) 



The state annihilated by all annihilation operators is an 
approximate ground state 



a a \0,t)[ ] = 0. 



(71) 



The vacuum expectation value of the Hamiltonian for 
fluctuations is 



R 3 



[0] (0,t\H ftG \0,t) [0] = —J2 



2 

3A / 1 



R 2 



/3 



The number states are similarly constructed by 
K,i)[o] = -i=[alj \0,t) [Q] . 



(72) 



(73) 



The quantum state of each mode is a linear superposition 
of its number states. The quantum state of the fluctua- 
tions is a product of such quantum states 



*/,[o](0/)=n*«. 



(74) 



As emphasized before, the LN approach constructs the 
Hilbert space for the classical background field and fluc- 
tuations on which various initial conditions can be rela- 
tively easily incorporated. This LN approach is expected 
to be a powerful and convenient method especially for 
studying the early Universe and the time-dependent pro- 
cess of critical phenomena. 



C. Comparison with the Gaussian Effective Potential 

In order to compare with the results from other meth- 
ods, let us consider the Minkowski spacetime 



R(t) = 1. 



The solutions to Eq. (|69|) are easily found 

1 



Vol 



-in a t 



The Q a satisfies the gap equation 

n 2 a = k 2 + ^ 2 (0 c ), 

where 

^) = m2 + 3A ^ + 3A (2^^i; 



(75) 



(76) 



(77) 



(78) 



is the renormalized mass. From the equations of the mo- 
tion, one is left with the quadratic Hamiltonian 



H 



a 

+ 3A 



3A0 2 )<p* a <p 



R 2 



(27T) 

3Ai? 3 1 



4 (2 



(79) 



First, we compute the effective potential at the zero 
temperature. The effective potential is obtained by tak- 
ing the expectation of the Hamiltonian (|5?]) except the 
term 7r 2 : 

m 2 A 3 A 

Kff.(^,M) = — 4>l + + Ji(M) - T Io(M), (80) 



where 



1 1 f d 3 k 1 



2 J (27r) 3 fV 
d 3 k 



(81) 



We have thus shown that the effective potential is equal 
to the Gaussian effective potential || . 

Next, we compute the effective potential at finite tem- 
perature. The effective potential equals to the Helmholtz 
free energy 



F 



HZ) = v ei 



(82) 



where Z is the partition function 
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Z = Tre 



-0H 



(83) 



To the first order of the coupling constant, the partition 
function is given by 



= Tre _/3 ^ G - A/3Tr[iJ}e _/3J?/ ' G 
The finite temperature effective potential is found 



(84) 



Kff.(0cM) 



A 



3A 



(Io 



T0\2 



(85) 



where 



1 



1 



1 



(2tt) 3 ^ fi a e^" - 1 
d 3 fc 1 1 



(2tt) 3 n(k) e^( k ) - 1' 



d 3 fc 
(2V) 3 



In 1 



,-/3H(k) 



(86) 



The finite temperature effective potential can also be 
rewritten as 



V e lX<p c ,fi) = V eS X<l>c,M) 



d 3 k 

f3j (2T) 3 



In 1 



-/3Q(k) 



(87) 



in terms of the rcnormalized mass 



M 



+ 3A(/ + IS + €)■ 



Hence the effective potential at finite temperature is 
equal to the Gaussian effective potential at finite temper- 
ature. We have shown that the LN approach to quantum 
field theory at the lowest order of the coupling constant 
gives the same results as those of the Gaussian effective 
potential approach. The real vantage of the LN approach 
is that one can compute readily the effective potential be- 
yond the Gaussian effective potential to any desired order 
as will be sketched in the next section. 



V. BEYOND THE GAUSSIAN APPROXIMATION 

We introduce an algebraic, nonperturbative method 
to quantum field theory. As a (0 + l)-dimensional toy 
model for a self-interacting quantum field, we shall con- 
sider again the anharmonic oscillator 



H 



1 



A 



(89) 



One of results of the NL approach is that the anharmonic 
oscillator has an algebraic structure which has not been 
much noticed in the literature. In the Fock space we can 
represent the Hamiltonian in the form 



H = QgA^A + E v 



(90) 



The operators A and do not satisfy the standard com- 
mutation relation. They rather satisfy a generalized de- 
formed algebra of the form 



AA* = F(A*A). 



(91) 



The standard commutation relation is given by F(x) = 
x + 1, and the g-deformed algebra by F(x) = Q 2 x + 
1. In the Fock space representation F is in general a 
polynomial of x. 

It has been found that up to the first order of the 
coupling constant the anharmonic oscillator has the su- 
persymmetric form || 



H 



n 



(4: 



(92) 



where the frequency is determined from the gap equation 

3A 



[i] 



(93) 



The explicit form of A and A^ is found from the LN equa- 
tions. Remarkably these operators satisfy the g-deformed 
algebra 



where 



Q 



1 



3A 



(94) 



(95) 



The first improved vacuum state beyond the Gaussian 
approximation is defined by 



(96) 



The g-deformed algebra enables one to find the excited 
states 



1 



n \i] = 



A [i] 



|0)[ 



where 



}2n 



(97) 



(98) 



It is straightforward to obtain the energy eigenvalues 

E [lln = ^(ln} + [n+l}). (99) 

To any order of the coupling constant, we are able 
to find the Hamiltonian in the form (|S9|). For instance, 



8 



to the third order of the coupling constant, we find the 
Hamiltonian 



H = Ml G A] 3] A m +E {3] 



(100) 



where 



= nn ( 



3 /h\ 
16 \W. 



3 (hX\ 2 

128 \n^J 



25 (KX\ 3 



1024 VftjU 



(101) 



-^[3], vac. i s the vacuum energy up to the same order of 
the coupling constant. The generalized deformed algebra 
( |9l| ) is found to be 

3 Ti A 

A m A\ 3] = l + A\ 3] A m + IW A\ 3] A [3] 



/?iA\ 2 r 9 
Vf^J I 



32 



g(4A) 2 } 



/?iA\3 r 387 141 Jt - 
+ (f^) 1 512 + ^4] ^3] 



+ S( i f3 1 ^) 2 + S(4 ] %) 3 }. 



(102) 



Therefore, one obtains the better improved vacuum state 
beyond the Gaussian approximation from 



^[3]|0)[8]=0. 



(103) 



As in the harmonic oscillator case, the excited states are 
obtained by applying A* to the vacuum state: 



l*)[3] 



1 



nto^ 



(A 3 ]) fc |o> 



[3] = 



where 



A [ 3 y 



(104) 



(105) 



]^[3]|*)[3] = JVfc|*>[3]- 

The eigenstates are constructed to be orthonormal 

[3](k\l)[3]=S kl . (106) 

It is straightforward to see that the eigenvalues satisfy 
the following recursive relations 

N k = F(JV fc _i), N Q = 0. (107) 

Therefore, the anharmonic oscillator has the energy 



E [3]ik = Ml G N k + E\ 



[3] ,vac. • 



(108) 



Though not shown explicitly, the algebraic structure 
of the quartic anharmonic oscillator with q-deformed al- 
gebra at the first order of the coupling constant or with 
the generalized deformed algebra at the higher order of 
coupling constant is expected to play a pivotal role in 



nonperturbative quantum field theory. For instance the 
q-deformed algebra results in the partition function of the 
anharmonic oscillator valid even for the strong coupling 
limit. The extension of the q-deformed algebra to the 
self-interacting quantum field in Sec. IV will give rise to 
much more correctly the effective potential at the zero or 
finite temperature than the Gaussian effective potential. 



VI. SUMMARY AND OUTLOOK 

In this paper we have introduced the so-called 
Liuoville- Neumann (LN) approach to nonequilibrium 
quantum fields. It is a canonical method that unifies 
the functional Schrodinger equation with the LN equa- 
tion. The LN approach is nonperturbative in that at the 
lowest order of the coupling constant it gives the same 
results as the Gaussian effective potential. The LN ap- 
proach is universal since it is equally applicable to time- 
independent and time-dependent quantum systems, and 
to equilibrium and nonequilibrium systems. Though not 
shown in this paper, the LN approach can be unified with 
thermal field dynamics by doubling the degrees of free- 
dom. This is achieved by introducing a fictitious Hamil- 
tonian H and enlarging the Hilbert space constructed 
in this paper. The LN approach combined with thermal 
field dynamics is expected to give more physical intuition 
to our understanding of nonequilibrium quantum fields. 
The renormalization problem and the partition function 
up to higher orders of the coupling constant are not dealt 
with. Interestingly enough, fermion systems and gauge 
theories are challenging problems to the LN approach. 



ACKNOWLEDGMENTS 

The author would like to thank Dr. D. Bak, Dr. J.-Y. 
Ji, Prof. S. K. Kim, Prof. K.-S. Soh and Prof. J. H. Yee 
for some parts of this work and many useful discussions. 
This work was supported in parts by the Center for The- 
oretical Physics, Seoul National University and by the 
Non-Directed Research Fund, Korea Research Founda- 
tion, 1997. 



[1] A. Linde, Particle Physics and Inflationary Cosmology 
(Harwood Academic Pub., New York, 1990). 

[2] For review and references, see N. P. Landsman and Ch. 
G. van Weert, Phys. Rep. 145, 141 (1987). 

[3] H. Umezawa, Advanced Field Theory (AIP, New York, 
1993); references therein. 

[4] S. P. Kim, Class. Quantum Grav. 13, 1377 (1996); in 
Thermal Field Theories and Their Applications edited by 
Y. X. Gui, F. C. Khanna, and Z. B. Su (World Scientific, 



9 



Singapore, 1996), hep-th/9511082 



" Liouville-Neumann 



Approach to Nonperturbative Quantum Field Theory", 
to appear in Recent Developments in Nonperturbative 
Quantum Field Theory, the Proceedings of AP CTP- 
ICTP Joint International Conference, Seoul 1997, jiep 



th/9706052| S. P. Kim, J.-Y. Ji, H.-S. Shin, and K.-S. 
Soh, Phys. Rev D 56, 3756(1997); K. H. Cho, J.-Y. Ji, S. 
P. Kim, C. H. Lee, and J. Y. Ryu, ibid. 56, 4916 (1997). 
[5] D. Bak, S. P. Kim, S. K. Kim, K.-S. Soh, and J. H. Yee, 
"Perturbation Method beyond the Variational Gaussian 
Approximation: The Liouville-Neumann Approach", 



hep-th/9808001 ; "Factorization and q- Deformed Algeb ra 
of Quantum Anharmonic Oscillator" , |iep-th/9808144 . 
[6] K. Freese, C. T. Hill, and M. Mueller, Nucl. Phys. B255, 
693 (1985); 0. Eboli, R. Jackiw, and S.-Y. Pi, Phys. Rev. 
D 37, 3557 (1998); J. Guven, B. Lieberman, and C. T. 
Hill, ibid. 39, 438 (1989); O. Eboli, S.-Y. Pi, and M. 
Samiulla, Ann. Phys. 193, 102 (1989). 
[7] H. R. Lewis, Jr. and W. B. Riesenfeld, J. Math. Phys. 

(N.Y.) 10, 1458 (1969). 
[8] S.-J. Chang, Phys. Rev. D 12, 1071 (1975); P. M. Steven- 
son, ibid. 30, 1712 (1984); ibid. 32, 1389 (1985); P. Cea 
and L. Tedesco, ibid. 55, 4967 (1997). 
[9] S.-Y. Pi and M. Samiullah, Phys. Rev. D 36, 3128 (1987). 
[10] D. Bak, S. P. Kim, S. K. Ki m, K.-S. Soh, a nd J. H. Yee, 
Phys. Rev. D 58 (in press), |gr-qc/98050"Tl| 



10 



